Invariant quantities in shear flow 
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The dynamics of systems out of thermal equihbrium is usually treated on a case-by-case basis 
without knowledge of fundamental and universal principles. We address this problem for a class of 
driven steady states, namely those mechanically driven at the boundaries such as complex fluids 
under shear. From a nonequilibrium counterpart to detailed balance (NCDB) we derive a remarkably 
simple set of invariant quantities which remain unchanged when the system is driven. These new 
nonequilibrium relations are both exact and valid arbitrarily far from equilibrium. Furthermore, 
they enable the systematic calculation of transition rates in driven systems with state-spaces of 
arbitrary connectivity. 
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Despite the abundance of nonequilibrium states of 
matter in nature, a fundamental understanding of the 
variety of phenomena exhibited in such states is still far 
from complete. There have been numerous recent theo- 
retical advances in non-equilibrium statistical mechanics, 
with a large body of literature devoted to establishing el- 
egant methods for calculating the macroscopic behaviour 
of stochastic systems, given their microscopic transition 
rates, and demonstrating the efhcacy of those methods 
by application to idealized models Instead, we ad- 

dress the prior question of how to choose the microscopic 
transition rates in the first place, when defining a non- 
equilibrium system. The simplest generalization of ther- 
mal equilibrium, namely nonequilibrium steady states, 
have many features in common with equilibrium steady 
states, and a mathematical description starting from first 
principles is possible in both cases. For a common class of 
systems — complex fluids in continuous shear flow — the 
nonequilibrium noise imposes certain rules on the rates, 
akin to (but different from) the equilibrium principle of 
detailed balance. In this letter we derive a remarkably 
simple set of invariant quantities valid in steady states of 
sheared fluids. On the basis of these new exact nonequi- 
librium relations we are able to formulate a systematic 
method to determine the driven microscopic dynamics in 
state spaces of arbitrary connectivity. 

We consider a probabilistic treatment of the system 
under consideration (a complex fluid) in which the dy- 
namical evolution is described in terms of a set of tran- 
sition rates {^ij} expressing the conditional probability 
per unit time to perform a transition between two mi- 
crostates i and j, where i,j G {l,...,rt}. We are inter- 
ested in the properties of the system when it has set- 
tled into a nonequilibrium stationary state under the in- 
fluence of continuous external driving at some distant 
boundaries. In the absence of driving, the set of tran- 



sition rates is constrained by the principle of detailed 
balance, stating that, for every pair of microstates i and 



(1) 



where Ei is the microstate's energy and the recip- 
rocal temperature parameter. Detailed balance results 
from four fundamental properties characterizing system 
and heat reservoir (i) ergodicity, (ii) microscopic re- 
versibility, (iii) stationarity and (iv) conservation of en- 
ergy. For a system with M transition rates there are 
M jl constraints in the form of Eq. ([1]) resulting from the 
statistics of the heat reservoir and ensuring equilibrium 
properties. 

Away from equilibrium, detailed balance is generally 
violated. In fact, the violation of Eq. ([1]) is often con- 
sidered to be the defining property of a nonequilibrium 
state. It was recently shown 0, [E, 0] that, for a partic- 
ular class of steady states including complex fiuids un- 
der shear, a nonequilibrium counterpart to equilibrium 
detailed balance (NCDB) can be derived. The crucial 
observation here is that, under continuous shear, a vol- 
ume element in the bulk of the fiuid has the same equa- 
tions of motion as at equilibrium; its Hamiltonian is not 
altered by the driving, but the stochastic infiuence of 
the surrounding fluid is no longer that of an equilibrium 
heat bath — the reservoir is itself under shear. As a re- 
sult the same properties (i) — (iv) apply as in equilibrium, 
amended only by an additional conserved quantity: the 
net average shear rate of all the fluid elements. Follow- 
ing either information-theoretic [1, 0] or Gibbsian Q ar- 
guments this small modiflcation to the ensemble results 
in an exact one-to-one mapping between the transition 
rates at equilibrium and those in the sheared steady state. 
That ono-to-one mapping implies a set of constraints on 
the transition rates in the driven ensemble. 

For continuous-time dynamics, the one-to-one map- 
ping of NCDB is expressed as [5] 



'Present affiliation: The Rockefeller University, 1230 York Avenue, 
New York, NY 10065, USA. E-mail: abauleOrockefeller.edu. 



(2) 



2 



Here, {r2y(i/)} denotes the set of transition rates in the 
sheared ensemble, parametrized by v, a Lagrange multi- 
pHer characterizing the driving strength of the reservoir. 
The LOij are the associated equihbrium transition rates 
that satisfy detailed balance, Eq. Under the influ- 
ence of shear forces applied at the boundaries of the whole 
ensemble, the rates are thus enhanced or attenuated with 
respect to equilibrium. The factor e^'^^J' exceeds unity 
if the transition i ^ j involves a conformational change 
that increments the shear strain by a positive amount 
Axij. This factor simply boosts every transition in the 
forward direction irrespective of the state space struc- 
ture. By itself, it would represent a simple mean-field 
expression for the driven transition rates. Important non- 
mean-field information about the global properties of the 
state space is contained in the quantities Aqji{i') which 
can be formally defined Q as 



Aqji{v) = lim In 



Eoo 
r— - 



G,(r,T)e- 



Er=-ooG'»(r,T) 



(3) 



The {Gi{r,T)} denote the equilibrium Green's functions, 
i.e., the probability that the equilibrium system, given 
it is in state i, exhibits shear strain r in time r. The 
quantities Aqji thus provide a measure for the increase 
(or decrease) in probability that the system exhibits the 
required shear rate (specified by v) at equilibrium if it 
performs the transition i ^ j. Therefore the mapping 
Eq. ^ expresses the fact that the likelihood of a tran- 
sition not only depends on the immediate flux contri- 
bution, but also on the new microstate's prospect for 
future flux. It is in this sense that NCDB goes sig- 
nificantly beyond simple mean-field expressions for the 
driven transition rates. However, the crucial ingredient 
in the mapping, namely the set {Aqji}, requires knowl- 
edge of all equilibrium Green's functions of the system. 
This requirement has hitherto impaired the practical ap- 
plicability of the theory. Below we demonstrate that the 
{Aqji} are in fact intrinsically related to the structure of 
the state space and can be determined following a sys- 
tematic set of network rules. We furthermore eliminate 
all Aqji values to derive some remarkably simple, exper- 
imentally testable relations between the transition rates. 

To this end we first introduce a graphical representa- 
tion of the state space following [3]. Here, vertices are 
assigned to the different states i of the system and edges 
to the possible transitions. If a transition is physically 
allowed to take place, i.e. toij > 0, then equilibrium de- 
tailed balance demands that the reverse rate ujji is also 
non-zero. Only connected graphs are considered, since 
they cover most physically relevant situations. In order 
to guarantee that the system can exhibit a macroscopic 
steady state, we require the state space to have a peri- 
odic structure along a coordinate x that measures shear 
strain. 

We define the basic graph as the graph corresponding 
to the non-periodic set of n distinct vertices. Its set of 
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FIG. 1: Example of a basic graph for a three-state system. 
The dotted line denotes an exterior edge connecting states in 
different periods. Here, states 4 and are of the same types 
as 1 and 3 respectively. E is energy and x is shear strain. 



edges will be denoted interior edges in order to distin- 
guish them from exterior edges connecting vertices of the 
basic graph with vertices of the next or preceding period. 
See Fig. ([T]) for a depiction of such a basic graph. The 
total number of transition rates in the system is '^i' 
where di is the degree (or connectivity) of the ith vertex 
of the basic graph including exterior edges. The minimal 
number of transition rates in an n-state driven system is 
2n, corresponding to a graph in the form of a simple con- 
nected path. For this class of state spaces the problem 
of finding the driven transition rates has a particularly 
straightforward solution due to the presence of invariant 
quantities in the steady state, derived as follows. 

Consider a system in microstate i, where it remains for 
a random duration t until making a transition to one of 
the di connected states. For a continuous-time Markov 
chain this waiting time has the exponential distribution 
hi(t) = (7^6"'^'*, where the total exit rate is defined as 



} 



The probability that the particle jumps to 



site j is then P^j — ujij/ai. The Green's function for state 
i can be determined as follows. From state i the system 
can only perform a transition to one of the connected 
states j; its subsequent displacement is then determined 
by the Green's function of state j. Gi is therefore related 
to the di Green's functions of the neighbouring sites by 



G,(r,r) = a 



/ di /i, (t - t ) V P„- Gj (r - Axj, , t) , (4) 



where the sum is taken over the set of states {j} con- 
nected with i. For normalized Green's functions, the co- 
efficient is Ci ~ 1/(1 — e^"'''^). We next transform the 
Green's functions to the corresponding quantities 



miii^,T 



In 



E 



G.(r,T) 



(5) 



whereby the difference rrij — is ultimately related to 
the Aqji in the long time limit (cf. Eq. ([3])). Multiplying 
Eq. Q by e'^^ and summing over the displacement r from 
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state i yields, after a shift in the summation variable, 
Substituting for Q, Pij and hi{t), and rearranging yields 



nii (i/,r)+o 



dte'^ 



(7) 



The integral can be removed by taking the derivative 
with respect to r on both sides, eventually leading to 



—mi[v,T)+a,-—mi{iy,T)e ' 

OT OT 



(8) 



In order to ensure steady state properties we have to 
consider the limit t — > cx) where rrij — mi ^ ^qji and 



d 



lim ^TOi(i^, r). 



(9) 



i.e., in the long time limit the rate of change of mi(z^, t) 
converges to a state-independent function of v, the flux 
potential Q{v) In Eq. ([8]), the result is 



E 



(10) 



On the right-hand side we can identify the transition 
rates in the driven steady state according to Eq. ([2]) • We 
therefore obtain a fundamental relationship between the 
equilibrium transition rates, the corresponding rates in 
the driven steady state, and the flux potential Q{v): 



Q{v) = - ai, 



(11) 



where I]i(z^) = J2 {j} ^ij i'^) ■ ^1- d) states that, for 
every microstate i, the total exit rate in the driven en- 
semble differs from its equilibrium counterpart by a flux- 
dependent but microstate-independent constant. On the 
basis of this central result a number of important impli- 
cations of NCDB can be derived. It turns out that it 
is not necessary to explicitly calculate the Green's func- 
tions of the system in order to determine the Ag's and the 
driven transition rates. Rather, as we will see more ex- 
plicitly below, the quantities of the NCDB formalism are 
intrinsically related to the graph structure via Eq. (jlip . 

It is now straightforward to formulate two sets of in- 
variant quantities for the steady state. The first was 
found previously for continuous-time dynamics [5] and is 
a consequence of the antisymmetries Aqjt ~ — Ag^j and 
Axji — —Axij. From Eq. 12]) the 'product constraint' 



directly follows: 



(12) 



Secondly, Eq. (|lip implies the ^exit rate constraint' 



E 2 ^ J — ^ i ^ 3 ' 



(13) 



NCDB thus predicts that the product of forward and 
reverse transition rates and the difference of total exit 
rates for every pair of microstates are the same as in 
equilibrium and therefore invariant with respect to driv- 
ing. No near-equilibrium assumptions have been made 
in the derivation, so the above relations are both exact 
and valid arbitrarily far from equilibrium. 

As well as being important and elegant in their own 
right, Eqs. and (fT5)) can be used to find the rates 
^ij (Eq. (12])) without evaluating the RHS of Eq. ([3]), as 
follows. Considering the whole basic graph with n states, 
there are n equations in the form of Eq. pip . This set of 
equations is sufficient to determine the unknown Ag^^'s 
in Eq. ([2]). Due to the relationship Ag^j = — Aqji, every 
edge of the basic graph is associated with two transition 
rates (forward and backward transitions) that depend on 
just one Ag. The number of independent Ag's is further 
constrained by closed paths in the graph, i.e. paths that 
begin and end at the same type of vertex, since the sum of 
Ag's along such a path vanishes (a 'loop constraint') due 
to Eq. ([3]). The total number of independent Ag's in the 
basic graph is always n — 1, as is seem by considering first 
the simplest basic graph, namely all n states connected 
as a simple path without any loops. In this case the 
number of edges is trivially n. Since one loop constraint is 
generated by the periodicity, there are n — 1 independent 
Ag's. From this simply connected graph all graphs of 
higher degree are generated by adding new edges. But 
adding an edge generates a new Ag and at the same time 
a new loop constraint, so that the number of independent 
Ag's always remains n — 1. With this knowledge, we can 
formulate the following network rules for the calculation 
of the driven transition rates in networks of arbitrary 
connectivity: 

Edge rule. Every interior edge and every pair of ex- 
terior edges in the basic graph corresponds to two rates 
that depend on one Ag in Eq. 

Vertex rule. For every vertex in the basic graph the 
difference between the driven and equilibrium total exit 
rates equals the flux potential Q (Eq. ITT]) ) 

Loop rule. For every closed path of edges the sum of 
the Ag's along this path is zero. 

In this formulation there are in total n equations and 
n unknowns, namely one Q and n — 1 Ag's. The number 
of independent equations can always be further reduced 
by eliminating Q, such that one is essentially left with 
n — 1 equations for n — 1 unknown Ag's. The solution 
of this system of equations fully specifies all the driven 
rates in the system as well as the imposed shear rate J, 
as functions of the flux conjugate parameter v which is 
related to the flux via dQ/dv = J [5]. 

Alternatively, an even simpler, t/-independent solution 
can be found from Eqs. (|12p and (fT5]) alone. However, 
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the number of constraints is then not sufficient to de- 
termine all the driven rates for every network structure. 
On the one hand there are Y^^=i ^^j/^ product constraints 
and n — 1 exit rate constraints. On the other hand, for 
an arbitrary graph configuration, there are X^iLi '^i tran- 
sition rates. Therefore only for graphs with the topol- 
ogy of a simple connected path (where di = 2), we can 
completely determine the rates using only the invariant 
quantities without applying the network rules. In this 
case we have 2n transition rates and 2n — 1 constraints 
stemming from the exact relations. The transition rates 
are fully determined if additionally the relationship be- 
tween the transition rates and the current is provided. 
In this formulation the driven rates depend on J directly 
instead of being parametrized by v. The relationship be- 
tween current and rates is further elucidated elsewhere 
0, but is straightforward to calculate for simple models. 

As an illustrative example we have applied the network 
rules to the three state system depicted in Fig.[TJ At equi- 
librium the system is completely specified by the set of 
eight transition rates which pairwise satisfy the detailed 
balance condition Eq. ([1]) . Let us now imagine that shear 
is applied to the system at the distant boundaries with a 
constant shear rate J. How does the continuous shear af- 
fect the microscopic dynamics of the system? This infor- 
mation is provided in a precise, unbiased way by NCDB 
applying Eq. If we know the shear contribution of 
each individual transition at equilibrium, i.e., the set of 
{Axji}, then the non-local quantities {Agji} can be de- 
termined via the above network rules. The results for the 
driven transition rates in the forward direction (positive 
Ax) are shown in Fig. [2l For forward shear three of the 
transition rates are strongly enhanced as expected and 
as would also be predicted by a mean-field model simply 
boosting each forward rate. However, we observe that 
the transition 1 ^ 2 is attenuated for increased shear, 
indicating that the system disfavors the path via state 2 
that requires two transitions to acquire the shear incre- 
ment Axai. For large shear rates the system will thus 
predominantly choose the direct path 1 — > 3 even though 
VLi2 also carries forward flux. A similar observation is 
made for the transition 5123 : for backward driving this 
transition remains significant because it connects to the 
favorable direct path. 

The search for fundamental principles governing the 
behaviour of systems in nonequilibrium situations has 
long been an area of intensive research. For systems 
maintained in a driven stationary state under the infiu- 
ence of a nonequilibrium heat bath, such principles can 
be identified in the form of NCDB, providing exact con- 
straints on the driven transition rates arbitrarily far from 
equilibrium. In this letter this theory has been investi- 
gated for systems evolving in discrete state spaces. As 
our main result we have derived a simple relationship be- 
tween the equilibrium and driven total exit rates which 
directly leads to a new set of invariant quantities. These 
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FIG. 2: The four forward transition rates of the three-state 
system of Fig. [T] with Ei = 3, £2 = 0, and £3 = 2 (/3 set to 
unity). Parameter values; ulJ\o = urj, = 1, iou = 2, 1^32 = 1.5, 
E2 = 0, Aa:;43 = 0.5, Aa;32 = 1.8, Aa;2i = 1. 



invariant quantities are non-trivial predictions of NCDB 
for arbitrary state spaces and can be considered as new 
exact relations in nonequilibrium statistical mechanics. 

Even in a simple three-state model, the non-local cor- 
relations of NCDB become evident. Whereas mean-field 
theories would simply boost a particular transition in the 
forward flux direction, NCDB takes into account a state's 
propensity for achieving fiux. A forward transition can 
thus be attenuated if it connects to a jammed state or 
from where subsequent transitions carry low flux. This 
striking property of the theory might ultimately be able 
to describe the counter-intuitive phase behavior exhib- 
ited for example in real complex fluids under flow. 
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